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O ■ ABSTRACT 

We consider resonant absorption in a spectral line in the outflowing plasma 
\^ I within several tens of Schwarzschild radii from a compact object. We take into 

OO I account both Doppler and gravitational shifting effects and re-formulate the 

<^ ■ 

[~^ . theory of P-Cygni profiles in these new circumstances. It is found that a spec- 

^— I ' 

• ' tral line may have multiple absorption and emission components depending on 

f^ ' how far the region of interaction is from the compact object and what is the 

f^*) , distribution of velocity and opacity. Profiles of spectral lines produced near a 

^ . neutron star or a black hole can be strongly distorted by Doppler blue-, or red- 

j^ ' shifting, and gravitational red-shifting. These profiles may have both red- and 

H ' 

C^ ' blue-shifted absorption troughs. The result should be contrasted with clas- 

sical P-Cygni profiles which consist of red-shifted emission and blue-shifted 
absorption features. 

We suggest this property of line profiles to have complicated narrow ab- 
sorption and emission components in the presence of strong gravity may help 
to study spectroscopically the innermost parts of an outflow. 
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1 INTRODUCTION 

It is widely accepted that radiation produced in the vicinity of compact objects can bear 
an imprint of the strong gravitational field, of either a stellar mass black hole (BH) or a 
neutron star (NS) in Galactic X-ray binary systems, or of a super-massive BH in case of 
active galactic nuclei (AGN). Thus the problem of direct detection of line emission from 
regions close to compact objects is of great importance. One of the most important piece of 
evidence for the presence of an accretion disk, and also the most extensively theoretically 
investigated subject in this field, is a broadened, red-shifted spectral feature attributed to 
a fluorescent Ka emission line observed in the X-ray spectra of some Seyfert 1 galaxies. It 
is believed to be formed due to the reprocessing of X-Rays emi t ted b y corona in the inner 



Fabian et al. 



(I2OO0I)). If this is the case. 



parts of a relatively cold accretion disk (see e.g. 
then the line forming region translates to a broad region in velocity and frequency space 
making a line profile severely skewed by relativistic effects of the accreting plasma. This 
strong distortion of the spectral line makes detailed diagnostics of physical conditions in the 
plasma using this line difficult. 

The grating spectrographs on the X-ray telescopes Chandra and XMM-Newton provide 
unprecedented spectral resolution up to ~ 10 keV. Thus, by means of Chandra High Energy 
Grating observations , a narrow "core" of the Fe Ka line was detected in many type I AGNs 
(JYaqoob et al.l (120031 )). In order to explain this spectral feature in these sources it is necessary 
to invoke both reflection from an accretion disk and also emission from distant regions (such 
as an obscuring torus). 

Observations of several quasars and Seyfert galaxies reveal narrow X-ray absorption 
lines indicating matter outflowing, falling onto, and orbiting near the BH. These features are 
interpreted as Fe absorption lines, blue-, and red-shifted relative to the rest-frame frequency. 

In some cases, observed spectral features are interpreted as being gravitationally red- 
shifted absorption lines. Tentative evidence has been found for an absorption feature which 
presumably originates from resonant scattering of iron within an accr etion flow in the stron g 



gravitational fiel d of the BH i n the Seyfert 1 galaxy NGC 3516 
In the paper by 



Nandra et al 



(|l999|)) 



Yaqoob et al.l (l2005l ) it was suggested that the detection of the narrow 
red-shifted absorption line superimposed on the red wing of a broad Fe Ka line in the 
quasar E1821+643 is due to resonance absorption by F e XXV or Fe XXVI , gravitationally 
red-shifted from region within 10-20 r^ from the BH. iReeves et al.l (120051 ) interpreted an 
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absorption features from the quasar PG 1211+143 as being produced by ionized iron and 
either being Doppler red-shifted due to the plasma falhng onto the black hole, or due to 



gravitational red-shifting, or both. 



Matt et al 



(I2OO5I ) also reported on the possible detection 



of a transient absorption line from the ionized Fe from the quasar Q0056-363. If the observed 
red-shift is only gravitational t hen it has been argue d the absorbing gas lies within ~ 5 r^ 



from the BH. In the paper by 



Dadina et al. 



(I2OO5I ) the red- and blue-shifted absorption 



lines are found in the X-ray spectrum of the Seyfert 1 galaxy Mrk 509. The data suggest 
these features are formed in a flow moving with velocity, v ~ 0.15 — 0.2 c, or/and being 
gravitationally red-shifted. 

Much smaller astrophysical objects which may potentially demonstrate gravitationally 
red-shifted lines are neutron stars. It is widely believed that X-ray bursts are produced as 
neutron stars are accreting plasma from their massive companions in close binary systems. 
It was also suggested that measurement of the gravitationally red-shifted absorption lines 
during an X-ray burst can yield the mass-to-radius relation for the neutron star and thus 



Cottam et al. 



provid es the required constraint to the equation of state of the NS interior. 
(I2OO2I ) reported on the observation of several absorption features in the X-ray burst spectra 
of the low-mass X-ray binary EXO0748-676. Prominent absorption features were attributed 
to the red-shifted lines of Fe XXVI, Fe XXV, and O VIII. It was speculated that both 
gravitational red-shift and possibly a slow outflow are responsible for the observed features. 

From the above examples one may conclude that there exists the possibility to observe 
gravitationally red-shifted narrow spectral features from diverse classes of objects. 

With an exception of a purely photospheric line, it is natural to assume that plasma 
which produces line emission is moving with respect to the source of radiation (photosphere 
or accretion disk). 

The line emitting plasma may be in the form of a diffuse wind and may or may not 
have a clumpy structure. In both cases one may hope to acquire important information from 
analyzing the proflles of such lines, and to learn about the strength of the gravitational fleld, 
distance, velocity fleld etc. 

Many normal, luminous stars possess strong mass loss. In many cases it is possible to 
derive information about the wind by analyzing the P-Cygni proflles formed in these winds. 
These peculiar proflles have long been known. The flrst paper w hich contained an e xplanation 



of P-Cygni proflles on the basis of the Doppler effect was that of 



Bealsl ( 11929. 193lh , who gave 



a basic explanation of the observed line proflles from novae and Wolf-Rayet stars. Since then. 
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these profiles have be e n stud ied both observationally and theoretically in numerous papers. 



For example, 



Morton! (119671 ). extensively investigated the P-Cygni ultraviolet resonant line 



profiles from winds of early supergiants. 



A theoretical breakthrough was made in the paper by ISobolevi (jl960l ) , in which he rec- 
ognized that if there exists a velocity gradient along the line of sight, the radiation transfer 
problem becomes purely local. That is, distant regions in the fiow can no longer exchange 
information using photons with frequencies within the width of the line. The radiation which 
is seen at a certain frequency within a line profile by a distant observer comes from a sur- 
face of equal line-of-sight velocity, (see Sect. 3.1). Following this, significant attention has 
been paid to computing and analyzing P-Cygni profiles and also to investigating their de- 
pendence on fundament al parameters of the fiow (see e.g. the atlas of P-Cygni profiles by 
Castor fc LamersI JiDTqI )) . 



The main goal in this paper is to calculate the line profile produced by a plasma which 
is moving in the strong gravitational field of a compact object. Thus, we want to obtain an 
answer to the question of how much the gravitational red-shifting can change the observed 
line profile in comparison with the standard case of P-Cygni line. Expecting that the grav- 
itational red-shifting may introduce some characteristic distortion to the P-Cygni profile, 
we wish to investigate what kind of the distortion is made and whether it can be used to 
deduce information about the wind and compact object system itself. 



The plan of the paper is as follows: in Sect. 2 basic assumptions are described about 
physical conditions and geometry of the line forming region. Here the objectives and limi- 
tations of the approach are defined; taking into account gravitational red-shifting, in Sect. 3 
we first calculate a spectral line optical depth. This is done in the Sobolev approximation 
which we adopt throughout the paper. In this section equal frequency surfaces (EFS) are 
re-defined to include the gravitational red-shifting effect and the shapes of such surfaces 
are calculated for several typical velocity profiles. The results obtained in this section are 
extensively used in the calculation of the radiation field within a spectral line in Sect. 4 or 
when we numerically calculate line profiles in Sect. 5. At the end of the paper we provide a 
discussion and summarize the results. 
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2 ASSUMPTIONS 

As mentioned above, the problem we address in this paper may be relevant to various 
astrophysical objects. Physical conditions in the wind driven from the accretion disk in 
AGN are different from those in the out-flowing plasma during an X-ray burst. At the same 
time we would like to elucidate some new features which may arise in a spectral line profile 
when the influence of the gravitational field is non- negligible. As a zero order model we adopt 
the simplest geometry and make further hypotheses about the velocity profile, temperature 
distribution, opacity etc. We wish to consider a minimum number of free parameters, as their 
overabundance will likely obscure rather than elucidate the origin of expected new features. 
We make the following assumptions about the geometry and kinematical properties of the 
absorbing plasma and the source of continuum photons: 

(i) The velocity distribution is spherically symmetric. This assumption allows a simpli- 
fication to the radiation transfer problem and also adds robustness to the results. In some 
cases a departure from spherical symmetry (e.g. wind from accretion disk) should be taken 
into account. For example, rotation may have some influence on the line profile. However, 
in radiationally or thermally driven accretion disk winds , the radial component of the ve- 
locity can quickly exceed the toroidal component, due to the conservation of the angular 



momentum. Notice, that this may not be the case in the MHD winds, in whi c h po 



oidal 



velocities may not be much larger than the toroidal ones e.g. iBlandford Sz Paynd (119821 ). On 
this ground we hypothesize that inclination and geometrical effects due to accretion disk 
may have more influence on the resultant line profile. 

(ii) The velocity profile is represented by a gradually increasing function of the radius. In 
the case of P-C ygni profiles from normal stars, decelerating winds w ere considered in several 



works ( see e.g. iKuan fc Kuhil (Il975l ). iMarti fc Noerdlingerl (119771 ) ). Some results of these 



papers are discussed further in the text. 

(iii) The only source of continuum photons is the spherical core. Accretion disk may 
influence the results in two ways: close to the BH, the disk provides strongly anisotropic 
radiation field; the disk may screen the emission from part of the wind. At the end of this 
paper we investigate the effect of screening by the disk which is viewed face on. 

(iv) No bending of the photon trajectories is taken into account. While this effect plays an 
important role in formation of lines from the inner parts of accretion discs (both Schwarzschild 
and Kerr) it does not play an important role in the "zero-order" model adopted in this pa- 
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per. The relative importance of bending can be inferred from the equation which describes 
the orbit of a photon in the gravitational field of a Schwarzschild BH: 

1 1 

— ~sinQ; + - — (3 + cos2q;), (1) 

where a is the polar angle of the photon's trajectory with the impact parameter b, Qq = b/vg, 

2GM 
Tg = — - — is the Schwarzschild radius, and x = r/b. In the current studies we are concerned 

with regions of the flow located approximately at radii > lOr^. Thus, taking x = oo, qq ~ 20 

we obtain the deflection angle, Aa ^ 5.7°. Taking into account bending of photon trajectories 

would introduce a significant complication in the solution the radiation transfer problem and 

would require frequency dependent Monte-Carlo simulations of the radiation transfer in the 

frame of the General Relativity (GR for short), which is beyond the scope of the current 

studies. 



3 OPTICAL DEPTH 

We associate an observer situated at infinity, 0°° with the laboratory frame ("lab"- frame for 
short). The fluid is moving with the velocity v. It is most convenient to measure absorption 
and emission coefficients in the co- moving frames, Ocom; coinciding instantaneously with 
the fluid at each point along the streamline. To obtain a relation between these two frames 
one needs to consider another local frame, Oioc - a Lorentzian frame which is at rest at a 
given point Sq and instantly coincides with the Ocom frame. Generally Oioc and 0°° are not 
the same because of the presence of the gravitational field. The opacity in the Oioc frame 
is obtained according to the transformation: xLc — Xcom '^/'^loo where Xcom ^^^ ^ ^^^ ^^^ 
absorption coefficient and the frequency of the radiation in the co-moving frame. Notice that 
we neglect all effects associated with the strong gravitational field except for the change of 
the radiation frequency. In particular, we assume that photons are propagating in straight 
lines, so that /iiab = /^loc, where /i is the cosine of an angle between the radius- vector and 
the direction of the ray. The frequency of a photon that is propagatin g in the background 



gravit ational field can be found from the relation: lyiocy/goo = const, (iLandau fc Lifshitz 
(1960)) 5 where f^oo is the corresponding component of the Schwarzschild metric tensor. In 



the "weak field limit" it follows that ^/goo ~ 1 + 0/c^ where 0(r) is the potential of the 
gravitational field. 



In this paper, we adopt the pseudo- Newtonian potential of Paczynski-Wiita (PW) (jPaczynski fc Wiita 

(|l98o|)): 
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GM 



rg — r 



(2) 



The PW potential mimics important features of exact general relativistic solutions for parti- 
cle trajectories near Schwarzschild black hole. This potential correctly reproduces the posi- 
tions of both the last stable circular orbit, located at Sr^ and the marginally bound circular 
his u seful property, to capture the essentials of GR effects, has been used in 



orbit at 2r^ 



Dorodnitsynl ([20031) to calculate the structure of the line-driven wi nd near compact obiect, 



and th e results are in a good agreement with the GR calculations of iDorodnitsyn fc Novikov 



(120051) 



The Lorentz transformations between Oioc and Ocom give: i) = 7z/ioc(l — /i/?), where 
l3 = v/c, 7 = (1 — /3'^)~^^'^. Consider a photon that is emitted at a point sq. At some other 
point s it's co-moving frequency is i>(s) = 'yi^iodso) J gooiso) / goo{s){l — fi{s)(3{s)). In order 
to allow for moderately high terminal velocities we will write all equations out to the second 



order, i.e. retain terms /3, /3^, and 0/c^. Taking into account that t'ioc('5o)Yfi'oo('So) = ^°°, 
where z/°° is the frequency measured by the observer 0°°, we obtain: 

oW = <'~(i-.W;3W-^ + ^). (3) 

The probability to emit a photon within a frequency range (z>, z/ + du), and within a 
range of solid angles in the co-moving frame {Q, Cl + dQ) is: dP^, = (4??)^^ \E'(z/ — z/q) dudQ, 
where uq is the frequency of the line in the co-moving frame, and the line profile function \l/ 
is assumed to obey the normalization condition: /q°° \1/(z> — z/q) du = 1. 

The optical depth between a point s and +oo along the photon's trajectory, which is a 
straight line in our approximation, is calculated adopting a transformation from space to 
frequency variable: 



~ / ,j r ''(^')..i 



t = I xLds'= ^^Xo,com* (^ - ^o) c?s' (4) 

Js Z/ioc 

i^(o<^) Z> I VP (z> - Z/q) ,~ 

Xoxom / 7~ \ "-^i 

D{s) Z/loc dl^\ 

[dsj^, 

where Xo com i^ the line-center opacity in the co-moving frame: 

I Tre^ Ni/gi-Nu/gu ,.. 

Xo,com = —{9f) ^^ , (5) 

where A^^, Ni and f7„, gi are respectively: populations, statistical weights of the corresponding 
levels of the line transition, / is the oscillator strength of the transition, Az/ = v^vthjc is 
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the Doppler line width, and Vth is the thermal velocity. Given the Sobolev approximation 
we assume that a photon, after being scattered in a line at a point sq, can further interact 
with matter only in the immediate vicinity of this point. Due to gradients of the velocity 
and gravitational potential, such a photon is quickly shifted out of the resonance with the 
line. We expand i) in the vicinity of Sq in Taylor series retaining two terms in (s — Sq): 

f(,) ^ f(,„) + gj (,_,„)+ Ijg (,_,„)^ (6) 

\ / So \ / so 

Usually in the Sobolev approximation only the first order term is left in (|6]), and in most 
cases this is sufficient. However, there are some pathological situations, namely a singularity 
in (j3]) in which the second order term in the equation ([2]) is required. 

Using relations (jlj) and ([6]) and taking into account that in the Sobolev approximation 
the opacity is considered to be constant throughout the resonant region, the expression (jl]), 
to the second order v'^/c^, reads: 

t ^ ro(/i, s) (1 + (3\1 + fi') - 2/X/5 - 4) C(^), (7) 

where 

^o(/i,s)~Xo,com/IQ(/^'^)l' (8) 

and 

Q(yU, s) = «! + ^2 As, (9) 



where 



-1 di>\ ^ _/,ooN-i/^^ 



and As = s' — s, where s' is a point located on the trajectory of the photon. In practice, 
we take As = const. This quantity is required in order to evaluate numerically the optical 
depth in those directions, yU at which dv/ds vanishes (see discussion further in this section). 
The second factor in (ITj), (^{u) reads: 

c(z>)=r ^{p-p^)dv. (11) 

For simplicity we assume that the line absorption coefficient is zero outside the frequency 
interval (z/q — Az//2, z/q + Az//2) i.e. in terms of Az/, spectral line has a half thickness 1/2. 
In equation ((Tj), it was assumed that (z/°°)"^ ~ {vq)^^{1 + /?^/2 — fif3 — 0/c^). 

We integrate the radiation transfer equation along characteristics, which in spherical 
symmetry are the rays of constant impact parameter, p. Differentiating in a direction of 



into account that 



ds 



dz , 
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ions of {r^n) or (p, -z). Taking 
z/ yjp^ + z'^1 and r = \/pi^ + z'^, 



s = zz, we consider all dependent variables being either functions of {r,^) or (p, -z). Taking 

d ( d\ 1 — fi"^ d d , 

+ /U7— , and /i 



r 



9/1 9r' 



we find: 



ly 



_^dv 
ds 



r 
+ /i 



/.Ml 



dlnp \ 
d\nr J 

d\n(3^ 

P 



(12) 



^P c^ dlnr dlnr 
Similarly, after some algebra we find: 



- [ly 



_^d^i> 
d^ 



fi^ f2fi(3 



+ /i 



r 
1 



fi 



\ r 



+ (3'{(3-2fi) 



(13) 



^-/?' +M(/9')'-/i/3" + /3/5"-^) 



r 



where the prime denotes differentiation with respect to r. 

In classical Sobolev theory, di) /ds appears to be an even function of /i. In our case a 

situation is possible where for certain values of /z Doppler shifting and gravitational shifting 

cancel each other, zeroing the right hand side of equation flT2l) . Note that even if the ^cf)' 

term in (TT2|) is taken away the right hand side of this equation can change sign depending 

on whether dv/dr is positive or negative if radius is increased. Thus for decelerating flows 

I-/.2 



it is possible that purely geometrical frequency shift. 



-V (due to the divergence of the 



spherically-symmetrical fluid flow), can compensate for the /i^f' term. 

Strictly speaking, for those directions of /i at which gravitational frequency shifting can- 
cels Doppler shifting, the first order (in (s — sq)) Sobolev approximation is formally not 
applicable because of the si ngularity foun d in the denominator of eq. (jHD- A similar prob- 



lem has been also found by iJefferyl (Il995l ) who considered the relativistic, time-dependent 
Sobolev approximation. In practice when numerically evaluating the integral (jl]) we adopt 
the following procedure: In most cases the first order Sobolev approximation works very 
well and only the term which is first order in (s — Sq) has been retained in ([8]). In our 
calculations, in most cases the singularity in (jlj) has not been detected. On the other hand, 
in those rare situations when it was found, we have taken the second order term into account 
in the denominator of ([7]). Doing so, we use a second order expansion ([6]), and additionally 
specify a constant value of As. It should be emphasized that our experiments in choosing 
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Figure 1. A photon is propagating in arbitrary direction FB from larger to smaller radii. A wind moves radially with the 
velocity V(r). Corresponding projections of the radial velocity on the direction of the ray at points so and s determine the 
relative Doppler shifting. For certain cases the latter can compensate for the gravitational shifting between these two points 
(see e.g. the discussion before and after eq. I|14p ). 



different values of As persuaded us that tfie resultant line profile is not influenced by this 
choice. 

To better understand the physical origin of the singularity, it is instructive to consider a 
ray formed by photons propagating from larger to smaller radii, i.e. from r(so) to r(s) in the 
direction FB as depicted in Fig. [H Assuming the velocity field to be spherically-symmetric 
and the velocity gradient positive, dv/dr > 0, and denoting /i(so) = cos 6*0, and fi{s) = cos 9, 
one can see that there are two cases to be considered: a) /i(so) < /^(s) < 0, so that 

IM^o)| > Hs)l (14) 

and b) /i(so) < , ^{s) > 0. Rewriting equation (l3l) and keeping only first order terms we 
obtain: 

6u/uo = -{(3{so)Hso)\TmHs)\) + ''^^'^' ~ ''^^""^ (15) 

where 5z/ ~ z/(s) — z/q, i^{so) — vo and the minus sign corresponds to the case a). If r(s) < r(so) 
and recalling that velocity is an increasing function of r, we obtain (3{sq) ^ I3{s). Since 
the last term on the right-hand side of flT^ is positive we see that in both cases it can 
compensate for the term which is due to the Doppler-shift. That is, in contrast to the case 
when no gravitational red-shifting is taken into account, here a photon emitted in the wind 
can again interact within the wind, thus providing a radiative coupling of distant points. 
Thus we arrive at the problem of multiple- valued equal frequency surfaces. 

The problem of multiple-branched equal velocity surfaces is known from modeling of 
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line profiles from decelerating atmospheres ( iKuan fc Kuhil (Il975l ). ISurdejl (Il977l )). These 



authors adopted a decelerating velocity law, v{r) = Vconst{Rph/ rY, where / > is a parameter 
of deceleration and -Rph is the radius of the photosphere. For such a velocity distribution a 
coupling of distant points due to the zero relative frequency shift is found. As a result of 
such coupling the radiation field at any of these points includes not only the contribution 
from the source of continuum radiation ("core") but also from distant parts of the flow. 

3.1 Surfaces of equal frequency 

It is convenient to introduce a non - dimensional frequency variable which measures a dis- 
placement of the frequency from that of the line center in terms of Doppler line width, 

Au: 

y = {u-uo)/Aiy. (16) 

Similarly, a co-moving version of this variable reads: y = {v — z/o)/Az/, and for the observer, 
0°^ we obtain y°° = (z/°° — z/o)/Az/ as a measure of the red/blue-shift within the line pro- 
file. According to the Sobolev approximation photons of certain frequency emitted in a 
spherically-symmetric, gradually accelerated atmosphere, may come only from the volume 
occupied by a thin shell of a thickness, l{p,z) ~ Vth/{dv/dr). In the limit, '^{y) = 5{y) this 
shell becomes a surface of constant line-of-sight velocity. This is a key idea behind the calcu- 
lation of the line profiles in the Sobolev approximation. The concept of surfaces of constant 
line-of-sight velocity, v cos 9 can be modified and expanded to account for gravitational red- 
shifting. In such a case they are be better referred to as equal frequency surfaces (hereafter 
EFS). Similarly to how the relation ([3]) was obtained, we write: 

y^-y = Av-\v^ -y)=^u + ^- wV(2C), (17) 

where u = v/v^h is the non-dimensional velocity, C, = c/vth, and $ = (f)/{cvth) is the non- 
dimensional potential. 

A photon of the emitted frequency uq (i.e. y = 0), emitted in the wind at a point (2:0, p), 
at infinity has, an observed frequency y°°. Then an equation for the resonant surfaces reads: 



y~ - ^^/ U [^^Zl + p2 j + $ ^^4 + p2 j _ _^2 ( ^,2 + p2 ) = 0. (18) 



This equation determines the locus of the equal frequency surface (i.e. Zq) as a function of 
p and y°°: The non-dimensional form of the PW potential, ([2]) reads: $(a;) = C/2(l — ^90)^ 
where x = r/R* is the non-dimensional radius, and R* is the radius of the spherical core 
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from which the wind is launched. A non-dimensional parameter, qq = R* /vg determines the 
relative importance of gravitational red-shifting (i.e. by equating go -^ oo one completely 
neglects the influence of the gravitational field on the energy of a photon). Equation (IT8l) 
may have multiple roots, ZQi depending on the values of y°° and p and the velocity law. For 
the range of parameters relevant to this work, it is possible to show that the last term in 
(TT8|) does not affect the overall properties of this equation. It is also true that for a given 
impact parameter, p, and for any reasonable velocity profile v{r) (i.e. smooth single-valued 
increasing function of r), the following cases should be considered separately: 

(i) Blue-shifting domain, y°° ^ 0: Equation (IT8|) may have zero or one root. The latter 
may be located only at positive z (which implies that if z ^ 0, the superposition of Doppler 
blue-shifting and gravitational red-shifting can result in the observable z/°°(?/°°) only for 
unique zq, p). 

(ii) Red-shifting domain, y°° < 0: The situation is more complicated: Equation flTSl) may 
have zero, two or three roots. At z ^ there is always one (gravitational red-shifting is 
stronger than Doppler blue shifting) or no root. At z < depending on the velocity law 
there can be a situation in which a superposition of gravitational and Doppler red-shifting at 
some radius ri (at which gravity is stronger but the velocity is small) can equalize the sum 
of Doppler and gravitational red-shifting at some larger radius, r2 > ri (i.e. where gravity 
is negligible but velocity is high). 

The shape of resonant surfaces vary depending on the velocity profile. We can understand 
their shapes by considering only a few characteristic cases for v{r). These cases are: 

(i) V = 0. This case makes sense if go ^ 0, i.e. photon is red-shifted by the gravitational 
field of central object. From flTTl) the maximum gravitational red-shift is y^i^ = (C/2) (1 — 
xgo)^^. For a given y°°, the EFS has the shape of a sphere of radius, x{y°°) = gQ^{l — 
C/(2t/°°)). 

(ii) V = const. An example would be the outer part of a stellar wind in which the flow is 
approaching terminal velocity moves at almost constant speed. Alternatively, considering a 
spectrum from a thin, moving shell one may approximate the distribution of velocity within 
such a shell as constant. In the case of negligible gravitational red-shifting, the condition 
Vz = fiv = const, where Vz is the velocity projected on the line of sight, determines the 
locus of the EFS. Thus, if f = const, the EFS has a conical shape. An intersection of this 
EFS with the p — z plane consists of two straight lines passing through the center and being 
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Figure 2. Equal frequency surfaces for the velocity law ||19|I . The parameters are: _R* = 25 rg and V°° = 0.01c {left) and 
R* = 10 rg and V°° = 0.1c (right). Notice that a surface of equal frequency can have several branches e.g. a ray with a given 
1^°° and p can intersect it up to three times (see text for details). The observer is located at 2 — ► oo. 



symmetrical about the z-axis. Even in the case oi Qq ^ and v = const there is a frequency 
shift along the line of sight originating from the divergence of the flow lines in spherically 
symmetrical geometry (see e.g. equation (TT^ ). 

(iii) V ~ const ■ r. Another name for this dependence is a "Hubble law", which describes 
explosive events, in which more rapidly moving particles at large radii are outrunning the 
slower moving ones at small radii. 

To understand qualitatively real winds, for example those which are quickly getting away 
from the potential well (for example as does a line-driven wind in O star) one may approx- 
imately consider their nonlinear velocity profile as being linear, t; ~ r in the inner part and 
V ~ const in the outer part of the flow. 
Fig. [2] shows equal frequency surfaces for the velocity law: 
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Figure 3. Equal frequency surfaces for velocity law: v{r) = V°°\jl , where V°° = 0.1c and R* = 15r 
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Figure 4. The same problem as in Fig. |3]but parameters of the flow are different: V°° = 0.01 c and R* = 15 ■» 



u{x) = U°^{x-l)/{xt-l), 



(19) 



where U°° = V^/vth, V°° is a terminal velocity, xt = Rt/R*, and Rt is the wind terminating 
radius, Rt = 100-R*. Thus, ii u{x) ^ U°°x/xt and go ^ oo (i.e. no gravitational red-shifting) 
one finds from flTSl) that EFS is determined by the condition z = Rt y°°/U°°. Notice that in 
Figl2] the equal red-shift surfaces for the blue-shifted photons are parameterized in terms of 
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?/max = U"^ (i.e the maximum blue-shift in the case of Qq -^ oo) . The shape of these surfaces 
resembles those obtained in the case of (70 ~^ oo. 

The most important difference from the case when (70 ^ 00 is the appearance of multi- 
valued surfaces at the single frequency, y°°. As was already mentioned, when y'^ < (red- 
shifting), equation flTSj) may have one root at positive z and two roots at z < 0. As 
an example, consider a case when y°° = 0.2 y^^^. There are two equal red-shift surfaces: 
an ellipsoid-like surface at small r and a plane-like at larger radii (Fig. [2|, left panel). At 
larger red-shift, y°° = 0.3?/min the inner surface shrinks to smaller r, where the gravitational 
potential is stronger, while the outer surface shifts to larger r, to the domain where the 
velocity and the corresponding Doppler red-shift are larger. 

Terminal velocity 0.01 c is two small to describe a continuous flow from such a small radii 
as R* = 25 rg. However, physical conditions in the flow (most importantly, the ionization 
state in the gas) should allow for the formation of lines from heavy ions and at the same 
time have enough column density to shield these ions from being completely ionized. In 
other words, when conditions to form the line are just right only in a given part of a flow, 
and where the velocity may be approximated by the linear dependence. For example, if the 
opacity peaks at small radii within a continuous flow, an effective maximum velocity in this 
shell-like region is much smaller than V°°. The low terminal velocity case is only applicable 
if the absorption takes place only at the low-velocity base of the wind. 

In a general case, we expect that a transonic stellar wind should have a terminal velocity of 
the order of the escape velocity at the wind launching point. For example, for the parameters 
in this paper, the escape velocity is of the order 0.2 to 0.3c. Fig. [2] (right panel) shows such 
EFSsfor V°° = 0.1c. 



(iv) V = V^Jl - — . (20) 

This is a typical velocity profile for the transonic wind that is launched deep in a potential 
well. Such a wind is characterized by a steep transonic region, containing a critical point or 
multiple critical points (as in the case of line-driven winds), and an extended plateau where 
the velocity approaches V^ . This type of a velocity law is usually considered to approximate 
gradually accelerated atmospheres of hot stars. 

The characteristic behavior of the EFSs in this case can be understood from the previously 
considered cases: v ~ const and f ~ r. Fig. [3] shows EFSs for the parameters V°° = 0.1c 
and -R* = 15 r^. One can see that the shape of the EFS is mostly affected by the outer part 
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of the flow, i.e. where v ~ V°° , rather then by the inner "acceleration" part. The latter is 
more important for the EFSs obtained for V^ = 0.01c (i.e. as in Fig. HI). 

Summarizing, we note that gravitational red-shifting adds new ingredients to the classical 
topic of surfaces of equal frequencies: 

(a) Equal frequency surfaces for the red-shifted part of the spectrum {y°° < 0) may be 
located both in front and behind the compact object. Equal red-shift surface can break 
suddenly as the impact parameter approaches some critical value. 

(b) For y°° < the EFS may have a closed shape (ellipsoid-like) form. A ray with impact 
parameter, p and frequency, y°° may cross resonant surfaces several times depending on the 
velocity law and on the relative importance of gravity, i.e. proximity to the compact object. 

Physical conditions (i.e. ionization balance, radiation field) which affect the formation 
of a particular portion of the line profile are very different on different branches of EFS. 
For example, in the simple case of i; ~ r the absorption part is affected by the EFS which 
resides close to the photosphere, while the emission part is created by the EFS that may be 
located quite far from the photosphere (because of the large surface area of such a surface). 
The same may hold true for some other velocity profiles (see section 5). Having said this we 
proceed to the next section and prepare final formulas for the calculation of line profiles. 



4 RADIATION FIELD IN A LINE 

The mean radiation field J at any point on the Equal Frequency Surface consists of photons 
emitted by atoms in the immediate vicinity of this point: a local contribution, Jioc, and a 
contribution from distant points, Jdist- The latter consist of photons emitted by the core, 
Jcore; ^iid photous whlch arrive from other branches of the EFS. We assume the intensity 
Ic emitted by the core is constant over the line frequency interval (z/q — Az//2, z/q + Az//2), 
with no limb darkening taken into account. 

In the presence of resonant surface Sj the calculation of the intensity is s traightforward. 



he n o n-local Sobolev approximat ion was mostly developed in the papers by (JGrachev fc Grinin 



fll975f ): iRvbicki fc Hummed (11978f )l 

After being emitted either by the core or re-scattered by another EFS the intensity 
remains constant along the ray s until it enc ounters the EFS at S j , whe re it changes discon- 



tinuously according to the relation (see i.g. iRybicki fc Hummeii (119781 )): 
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Figure 5. Illustration to the geometry of equal frequency surfaces (see text for details). Not to scale. 
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(21) 



The situation is illustrated in Fig. [51 In this figure S^ refers to the surface of equal frequency, 

and " — " or " +" sign denotes whether the corresponding part of the EFS is located at 2 < 

or at 2; > respectively. As was already established, the EFS may have a closed or open 

shape. Thus Fig. [5] represents a situation when gravitational red-shifting is important and 

there are two branches of the EFS: a closed, ellipsoid-like, and a plane-like at larger r 

(smaller z), c.f. e.g. Fig. [2l The second possibility is that the two branches of the EFS are 

both open surfaces. (To illustrate this point we draw such surface, S"^ with a dashed line). 

Other notation has the following meaning: SJ~ is a negative part of that branch which has 

a closed form and S2" is a negative and open part of the multiple-branch surface. 

The mean radiation field at the point 1 is contributed by photons arriving from the core 

and from surfaces S^ and Sg". Photons emitted by S^ if not occulted by the core, may be 

attenuated by the line scattering only locally to the point 1: J ~ / Ii{Vi)l3} dVl, where /i 

Jvt' 

is the radiation intensity coming from points on Sj" seen from point 1 within the solid angle 
VL' occupied by the unobscured part of Sf . One should note that the resonance surface 
which is seen from point 1, for example, Sj" is not T^i but a different EFS. In a general case 
it should be calculated separately for each point where J is calculated. The same remains 
true about all other possible branches of EFSs which can contribute to the mean radiation 
field at a given point. Generally speaking Fig. [2] only shows that screening is possible. The 
goal of taking all these geometrical effects into account is beyond the scope of this paper. 
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However, when calculating the radiation field at point 1 one may want to take into account 
contributions only from points 1' and 1", i.e. only from those points which are situated on the 
ray {p,y°°)- Their positions remain unchanged and this fact greatly simplifies calculations. 
This is the first simplification in our description of the radiation field. 

The second simplification arrives from the ignorance of possible mutual screenings of the 
EFSs between each other and between a core (see discussion below). The probability for a 
photon to penetrate to point 1 without being scattered in its immediate vicinity can be cast 
in the form: Pl^^ = (1 — e^'^^^ ))/r(S+). The intensity can be determined from the following 
relation: /i(l') = S{1'){1 — e^^^^'^), where S{1') is the source function at the point 1' on T.^. 
Analogously a contribution from 1" can be constructed. 

Radiative coupling between different branches of EFSs poses a difficult problem to treat 
self-consistently. To overcome this difficulty different authors adopt differe nt approximations. 



When considering radiative coupling between fiducial points s and s' 



Grachev fc Grinin 



( I1975I ) have ignored the variation of the sou rce function, assumi ng that S{s') = S{s). An- 



other approximate approach was adopted by iKuan fc Kuhil (119751 ) who ignored the radiative 



c oupling between s and s'. In the literature this is known as the 'disconnected approximation' 



Marti fc Noerdlinger 



(119771 )). Giving an extensive analysis of the problem of interconnec- 



tion of different EPS branches for decelerating fiow, these authors arrive to the conclusion 
that i) taking into account radiative coupling can be of importance for the calcula t ion o f 
the radiation force; ii) fortunately, the disconnected approximation of iKuan &: Kuhil (119751 ) 
which completely neglects such coupling, works well and gives generally correct results for 
the source function and resultant profiles. 

In this paper the appearance of additional branches of the EPS is treated in the spirit of 
the "disconnected approximation". Por example, when calculating mean intensity at point 
1" of PiglSl we neglect any contribution from those resonant points which could be located 
on the line from the core to this point. These could have important consequences because 
of the possible screening effect. In turn this would affect the emission part of the spectrum. 
To simplify the treatment we also assume that the core serves as the only source of photons 
which contributes to J at a given point. Por simplicity we assume that the core is radiating 
with constant intensity Ic- The intensity of the radiation at r ^ 00 at a given impact 
parameter p and with a given frequency y°° is 
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+ 5(S2-)(1 - e-(S2-))e-W^+)+-(^r)), 

(22) 

where r(p, y°°, E"*") is the optical depth on the corresponding branch of the EFS. The notation 
here is the same as in FiglS) S+ corresponds to the branch of the EFS which is located at 
z > 0. Accordingly, E,~ is located at z < and the subscript 1 corresponds to that branch 
which lies closer to the observer. 

Relativistic effects may be important when calculating the escape/penetration probabili- 
ties and the source function. However, the Sobolev approximation makes all radiation trans- 
fer purely local. How quickly a photon is getting out of the resonance is completely controlled 
and incorporated in to the optical depth. The major difference with the non-relativistic case 
is from the aberration effects. That is, the radiation field emitted isotropically by the core 
is not seen isotropic in the frame of the fluid. In the fluid frame the emission is assumed 
isotropic. Thus, in this frame, the expression of the source function, S will have the same 
look (in terms of penetration/escape probabilities) as in non-relativistic case. In the rela- 
tivistic case, the derivation o f this probabilities and the of the source function is given by 



Hutsemekers fc Surdejl (119951 ). 



Having already cal c ulated the optical depth, (expression ([8])) we use the arguments of 



Hutsemekers fc Surdejl (119951 ) to derive the probability of a photon emitted in a line transi- 
tion to escape in any direction from a given point in the envelope. To the order f ^/c^ it is 
expressed in the form: 

/Sesc = ^ £ A'oc (l + /?'(3/i' - 1) + 2/i/?) dn, (23) 

where /Jj^^ denotes the escape probability in the direction /i: 
_ (l-exp(-ro(/i,g)) 

To{s,fi) 

where tq is given by ([8]). The penetration probability in the frame of Sobolev approximation 
reduces to 

P^en = l f' PLd^^. (25) 
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where fic = cos 9c, so that 9c = arccosWl 



X'' 



is the maximum angle at which the core is 



seen from the point x. 

Weighting the intensity by the probabihty of a photon to penetrate to a given point we 
obtain the mean intensity, J: 

with Pi^c being defined by (1241) and ly is a dilution factor: 
1 .. 



(26) 



W 



s 



The source function in the case of pure line scattering, in the lab frame reads: 
J 



Pc. 



l+/3'(3/i' 



2/i/3 , 



(27) 



(28) 



where /3esc should be found from fl23|) and terms of the order v'^/c^ were retained in fl28l) . 



In the lab frame the source function (128!) is clearl y anisotropic , taking 



the direction of motion. For relativistic formula see 



non-relativistic, e.g. lMihalad (1l978l ). ICastorl (llQTOl ). A considerable simplification is obtained 



Hutsemekers &: Surdej 



l igher values in 
JlQQsk and for 



in the case of a linear velocity law v{r) r^ r, when from ( l28l) and ( l26l) . in non-relativistic 
case one finds: 



S, 



lin 



A 



pen 



where 



A 



pen 



Pesc Pesc 

Note that formulas fl29|) are applicable only in the case of a linear velocity law and when 
$ = at ([8]) so that Pi^^ does not depend on /x. In the case of an arbitrary velocity law 
or/and if the gravitational red-shifting is taken into account the integral in fl26l) must be 
evaluated in its general form. 

Here we again emphasize our approximation in which we ignore the influence of the 
additional branches of the EFS on the mean intensity, J. That is we ignore both negative 
~ e""^ and positive ~ ^(l — e^"^) contributions which may arise from resonant points on 
the line to the point where J is calculated. 

After the intensity, I°°{y°°,p) has been calculated from (|22l) . the normalized flux that is 
registered by the observer at infinity equals to 

F{y^,p)/Fc = r I^{y'^,p)pdp, (30) 

^ u 

where Fc is the fiux emitted by the core. 

When calculating the red-shifted absorption features, it has been found that the line 
profile at some frequencies displays notable oscillations. These oscillations are unphysical and 
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are due to the spec ific behavior of the source f unction on S+ (see Fig. [5]). The same effect 



has been found by iMarti fc Noerdhngerl (119771 ) for a decelerating wind from a normal star. 



We find that the amplitude of these oscillations is reduced (very slowly in our calculations) 
when taking more points pi. This is in fact a numerical artifact which is related to the fact 
that the brightness of S+ strongly depends on p when p approaches its maximum value 
Pi = max(p on S^). Thus, prior to calculation of the integral fl30|) we calculate pi, -the 

rpi rPout 

intersection of S+ with z = plane. Then we split (l30l) into / I°°pdp+ / I°°pdp. 

Jo Jpi 

This procedure allows us to eliminate completely spurious jumps of S on the boundary of 
S+. It can be also useful in the calculations of decelerating winds from normal stars where 
there are EFSs with sharp boundaries. 

5 CALCULATION OF LINE PROFILES 

Before we proceed to closer examination of line profiles we need to make several assumptions 
about physical conditions existing in the fiow. It is assumed that a spherically-symmetric 
wind originates at the photosphere which emits radiation in continuum. The radius of the 
photosphere is R*. The radiation emitted by the core is resonantly absorbed (scattered) in 
a line of rest frequency z/q in the moving plasma. The relative importance of gravitational 
and Doppler red-, blue-shifting is controlled by parameters R* (i.e. 5^0 = R* /^g) and \^°°, 



respectively. Addi t ional 



y, we specify the distribution of the opacity. Following the recipe of 



Castor fc LamersI (119791 ). we parameterize the radial optical depth t^ = t{fi = 1) (c.f. eq. 



([7])), as a function of the velocity. The following dependencies are considered: 

r^{w)=To{k + l)\-—^, (31) 

yi Wc) 

and 



k 
W 



rr{w)=Uk + l)^^—-^^, (32) 

where w = v/V°° is the non-dimensional velocity, Wc = v{R*)/V°^ = 0.01 and /c is a free 
parameter. The parameter Tq is related to the total optical depth at the line center. 

We adopt the following velocity laws: the linear velocity law (IT9l) . and the "stellar" type 
velocity law: 

/ 1 \°^2 

w = w,+ il-w,)[l-—j , (33) 

where ai, 02 determine the slope and shape of the velocity profile. The velocity and opacity 
distributions for different pairs of ai and 02, and for k = 2 and k = —2, are shown in Fig. [61 
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Figure 6. Velocity (left) and opacity (right) laws. Velocity law I I33I I with Wc = 0.01. Opacity curves: solid line: equation I I31II : 
dashed: equation l|32l l. k=-2; dot-dashed: equation (1321 . k=2. Curves are marked by pairs of a\,a2 from II33I I. 

For a given (70 = R* /"^gi "we specify ?/°° and calculate v{r) from (TT9l) or (!33|) . and then 
calculate equal frequency surfaces from flTSl) . then calculate r^ from fl3T|) . fl32l) . and the source 
function from (128!) . and then the spectrum from (l30|l . 

Given the complex shape of the EFSs, i.e. depending which branch of the EPS is tracked 
when looking for the roots of the equation flTSj) , we switch between p, 2; or r , /x as independent 
variables. Additionally, in all calculations presented in this paper the occultation by the core 
is taken into account. Other parameters of the model are Vth = 5.7 ■ 10^ cm s~^, Rt = 100 R*. 
The results are shown in Fig. [7] - Fig. [TT] and we consider them in turn. 



5.1 



r law 



In Sect. 3 it has been established that in the case of linear velocity law (1191) . the EFS (the 
red-shifted part) may have two disconnected branches if the influence of strong gravity is 
taken into account and provided parameter go is small enough (c.f. Fig. [2]). This is easily 
understood from the following arguments: in a normal stellar wind go = R* /vg — > 00, and 
the absorption feature may be formed only in that part of the wind which is approaching 
the observer, and thus is observed as blue-shifted. In case of a static configuration {v = 0) 
the EFS has the shape of a sphere of radius x = \(/2goy°°\. A segment of this sphere which 
is projected on the core forms the red-shifted absorption line. Red-shifted photons could 
come both from the receding part of the wind which is behind the core, and also from that 
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Figure 7. Profiles for R* = 25rg, V°° = 0.01c (upper left), R* = 25rg, V°° = 0.2c (upper rigfit), R* = lOrg, V°° = 0.1c 
(lower left), R* = lOr^, V°° = 0.3 c (lower right), and the linear velocity law I I19I I. Curves are labelled by the total optical 
depth To- Notice the different vertical and horizontal scales in different panels. The left hand plots have V°° < Vesc, where 
Vesc is the escape velocity at the wind launching point. Vertical axes: flux. Horizontal axes: frequency shift, y°° . 
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Figure 8. The effect of the distribution of different absorption laws: t,- ~ (1 — w) (left) and Tr ~ w (right). Profiles for 
R* = lOrj, V°° = 0.3c and To = 4 and the linear velocity law I I19I I. Curves are labeled by the parameter k. Vertical axes: flux. 
Horizontal axes: frequency shift, y°° . 

branch of the EFS which is in front of it. In some cases, the re-emission from the receding 
gas may be faint since the source function is smaU: S ~ r^^, and we may expect to observe 
a red-shifted absorption feature. In such a case the resultant profile will be characterized 
both with blue- and red-shifted absorption features. 

Fig- © shows profiles for different terminal velocities and different launching radii. The 
strongest second absorption is observed for the lowest terminal velocity V°° = 0.01c (upper 
left). Consider it in more detail: the maximum width of the red-shifted absorption is set 
by — l/max and ^min; in the picture, these are: ymax = 52.2 and ymin = —104.4. We can see 
that the width of the second absorption component is of gravitational origin. Parameter V°° 
determines the width and strength of the blue-shifted emission component. 
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Profiles for larger terminal velocity V°° = (0.1 — 0.3) c and different Qq show that for 
larger V°° and for a given y°° < 0, the EPS is located closer to the core and scatters more 
radiation, smearing the red-shifted absorption line. A peculiar characteristic of these profiles 
is a red-shifted absorption line superimposed on the background of the red-shifted emission 
component. Por example, if 5^0 = 10, V°° = 0.1c, (lower left), the red-ward edge of the 
emission is set by — ymax = —522 and the red-ward edge of the red-shifted absorption by 
l/min = -261. 

In PiglH] is shown the effect of different distributions of the opacity. Por r^ ~ (1 — w)^ 
(left), resultant profiles are generally similar to those oi k = 2 and r^ ~ w^ (right). However 
in other cases the results are quite distinct for different k. 



5.2 More realistic v{r) 

Pirst we calculate line profiles for parameters of the flow relevant to tho se of a wind from 



a norm al star. These provide a good test of our results against those of 



Castor fc Lamers 



( 1l979h . Results are shown in PigM They are calculated for different values of Tq. The terminal 
velocity is V^ ^ 1897 km s~^. 

The parameter qq = 25 ■ 10'^, is taken large enough to provide that gravitational red- 
shifting is completely negligible in this case. Thus one expects to see no deviation of the 
resultant profiles from those of P -Cygni. The profiles show n in Pig. [9] are in good agreement 



with those presented in Pig. 7 of ICastor fc LamersI (119791 ) (note, however that these authors 



use v/V^ as a frequency variable). In this case, the wind approaches the terminal velocity 
much more quickly than in the previous two cases. Profiles for the velocity law ( l20l) are 
shown in Pig. iQ (right) for Tq = 4. The increasing influence of the gravitational red- 
shifting (decreasing oi go) results in shifting of the edge between the emission and absorption 
components to the left (to larger red-shifts). Por example, at Qq = 15, the position of the 
edge ?/°° c^ —155 is roughly determined by the E^ component of the EPS (c.f. Pig. [3] at 
y°° = 0.8,?/min = 140). An interpretation is that the red-shifted absorption almost exactly 
compensates for the red-shifted emission. Increasing go from 15 to 250, results in a shifting 
of this edge to the right, finally producing line profile of the type shown in Pig. [9] (left). 

As in the case of li near velocity law, different distributions of opacity can significantly 
change the results (e.g. ICastoii (Il970l ) ). The opacity is modified by varying the parameter k 
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in (13T!) . (!32|) . Additionally, parameters ai and 0:2 can also be varied if using a general form 
of the velocity law ( 133|) . 

A set of profiles for r ~ (1 — w)^ is shown in Fig. [TOl The profiles are calculated for 
R* = lOvg, V°° = 0.3 c, and Tq = 3. Each panel shows them for different pairs of ai, 0:2 for 
the particular value of k. Again for ai^2 7^ 1 a narrow absorption line is superimposed on 
the broad blue-shifted emission line. Notice, that in Fig. [101 the slowly accelerating wind 
(^1,2 = 1) shows more absorption, than in other cases of steeper accelerated winds. The 
center of gravity of the emission peak in this case is dominated by photons coming from 
EFSs which are located quite far away, behind the core (wind is accelerating to slowly). These 
surfaces reflect too few photons {S ~ 1/r^) and as a result, the gravitationally red-shifted 
absorption line dominates over the red-shifted emission. 

Fig. fm shows the results for the same set of parameters as before but now for r ^ w^ . We 
see that the results are considerably different from the previous case. Only for 0:1,2 = 1 are 
there profiles which look like classical P-Cygni profiles. In most cases, the narrow absorption 
is superimposed on the blue-shifted emission. For k = —3 and k = —2, profiles look similar 
to those obtained from the linear law (c.f. Fig. [7]). 



5.2.1 Shielding by the disk 

Our method, as presented in this paper, does not allow to consider anisotropic radiation 
field of an accretion disk or some additional attenuation of the emission. However, we can 
gain some insight by considering a wind that is viewed face on (i.e. perpendicular to the 
disk plane) and accounting for the blocking of those photons which are coming from behind 
such a disk. Since a significant part of the emission profile is formed by these photons, we 
want to know what happens if this emission is reduced. Thus we attenuate those photons 
by a factor of e^'^''^, where Tgh is the attenuation optical depth. The results are shown in 
Fig. [121 As may be expected from the previous results, suppression of the emission leads to 
absorption-emission-absorption profiles regardless of the velocity law. That is because the 
emitting branches of the EFS, i.e. Sf , Sg" in Fig. [5] are screened but the absorptive branch 
S+ is still there and absorbing radiation. 
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Figure 9. Profiles for v(r) = V°° (l - -Sl^. On the left the wind is launched at R* = 25 • 10^ r^ , having V°° = 1897 kms"!. 
Each curve on this panel is labeled by the total optical depth, Tq. On the right the wind has V°° = 0.1c and Tq = 4, and 
curves are labeled by the parameter go = R* /'''a- In ^H cases the optical depth i ~ (1 — w). 



6 DISCUSSION 

The strongest double absorption profile is found for the f ~ r velocity law. The superposition 
of the gravitational red-shift and Doppler shift easily produces two pronounced absorption 
features if the velocity law is not steep (as in the case of linear law) or the maximum velocity 
within the line forming region is small enough (as in Fig. [TJ upper left). More realistic 
velocity law (1551) also gives W-shape profiles for several different values of ai, 02. In case of 
a greater terminal velocity v = 0.1 — 0.3 c (i.e. of the order of the escape velocity at the base 
of the wind) in most cases the absorption feature is superimposed on the red emission wing 
of the line for both r^ ~ (1 — w)^ and r^ ~ w'' distributions of the opacity. For the latter 
distribution, the W-shape profile is most pronounced, sometimes having approximately equal 
absorption troughs in red and blue part of the spectrum (Fig JTT]) . Blocking of the emission 
(as in the case of an obscuring disk) also results in W-shape profile (see Fig JT^ . 

In order to produce an absorption feature, a certain resonant surface (or surfaces) must 
reside on the line of sight between the source of photons and an observer at infinity. In the 
case of the stellar (low gravity) wind and provided the velocity is gradually increasing, each 
EFS is single- valued. On the other hand, if the wind is decelerating, the EFS is not unique. 
If this is the case then photons emitted by one branch of the EFS can be scattered by the 
other and vice-versa. As was shown above, the gravitational field introduces several new 
aspects, changing both the shape and the locus of EFSs. The most important difference is 
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Figure 10. Profiles for R* = lOrg, V°° = 0.3c, Tq = 3. The velocity law II33I I and r '^ (1 — ui)*. Legend; curves are labelled 
by pairs of parameters ai, Q2- 

the appearance of the resonant surface in front of the core and which is located in the red- 
ward part of the frequency domain. This component of the absorbing surface intersects the 
hne of sight to the core and scatters any hne radiation that crosses it. While multi-branched 
EFSs are not new, the appearance of the surface of equal red-shift in front of the core for 
outwardly accelerated flow is not known in the theory of stellar winds. 

This may be made clear if we imagine a static and extended distribution of plasma 
around a compact object. In this imaginary exercise, there is only gravitational red-shifting 
which is important and the corresponding EFS is just a sphere that is concentric with the 
center of gravity. 

In a moving plasma this spherical EFS is transformed depending on the velocity law 
such that for certain red-shifts there is a component in front of the the core (c.f. Fig. ([21)). 
Absorption within this EFS may form an absorption feature in the red-shifted part of the 
line profile. 
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Figure 11. The parameter setup and labeling of the curves are the same as in Fig. |9] The opacity distribution: t ~ w^ . 



The emission component is formed in the same way as in stellar case, namely because 
of the large surface area of those EFSs which are not intersecting the line of site between 
the observer and the core. This red-shifted emission is superimposed on the absorption line. 
The net result depends on the distribution of the opacity, i.e. on the amount of scattered 
radiation. The blue-shifted part of the profile is formed because of the interaction of the 
blue continuum photons with resonant EFSs giving the absorption feature for the blue part 
of the profile. The value of V°° roughly sets the width of the red-shifted emission and blue- 
shifted absorption. Results show that the centroid energy of the blue-shifted absorption 
depends, on the distribution of the opacity, and can be seen at significantly smaller shifts. 
If the gravitational red-shift effect is significant but smaller than the pure Doppler effect, 
there exists some maximum red-shift at which the EPS has a component at 2; > 0. This 
component attenuates the core radiation as Ice~'^^^ ^ ( 122!) . and produces an absorption 
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Figure 12. The effect of shielding of the emission component. Curves are calculated for r ~ (1 — ui); left: Tgh = 1; right: 
f"sh = 10; labeling: total optical depth Tq. Other parameters: _R* = lOrg, y°° = 0.3c. 



feature superimposed on the emission component, as in Fig. [T0lfT2l . Here, again, the details 
depend sensitively on the distribution of the opacity. 

The intensity I{y°°,p) is composed of the radiation of the core, ~ IcG~^ and the contri- 
bution f S{x,y){l — e~'^i) dxdy, where the integral is carried along the resonant surface. 
In the case of pure line scattering: S ~ 1/r^, and the "brightness" of the emission line is 
determined by the surface area of the resonant surface, balanced by the 1/r^ divergence of 
the radiation flux. Thus, it is not enough to have opacity concentrated to smaller, or larger 
r. If at small r the absorption and gravitational field are strong the absorption feature may 
be strong but the red-shifted emission formed by scattered radiation is also strong. In the 
other extreme, at large r the emission may very faint because of the small flux. The shape 
of the EFS is also important. Without actual calculation of the EPS there is no way to tell 
what is stronger: the red-shifted emission or red-shifted absorption. 

The conclusion is: if absorption takes place in a spectral line within a wind provided that 
the plasma is moving radially and gradually accelerated in strong gravity, a profile with two 
absorption features is quite possible; in some cases there is a prominent emission component, 
and the profile is observed as W-shaped. 

The results are sensitive to the assumptions about the opacity and velocity; the geometry 
of the EFSs depends on the distribution of v{r) and on g^ = R* /vg] the amount of scattered 
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radiation depends on the distribution of the opacity. If it is concentrated towards the star 
(as in Fig. EI), the conditions favor the formation of a W-shape profile. However in each case 
different possibihties in choosing v{r) and Tr{r) should be considered. 

7 CONCLUSIONS 

Our goal in this paper has been to study shapes of spectral lines from plasma which is rapidly 
moving in the vicinity of a neutron star or a black hole. The latter case can account for both 
stellar (Galactic Black Hole Candidates) and supper-massive BHs (AGN). In the well studied 
case of winds from normal hot stars, a rapidly moving wind interacts with the continuum 
radiation of a star and produces a P-Cygni profile. From the theory of stellar winds it is 
known that in the case of an arbitrary spherically-symmetrical distribution of plasma which 
is moving with a gradually increasing velocity, an absorption line that is blue-shifted with 
respect to the emission line is observed. This is widely interpreted as a fingerprint of moving 
plasma in a variety of astrophysical situations. 

In this paper, shapes of line profiles for several velocity and opacity distributions were 
calculated taking into account gravitational red-shifting. 

Strong gravitational red-shifting helps the radiation to escape efficiently and to interact 
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(I2OO3I ) and IDorodnitsyn &: Novikovl (l2005l ) addressed a problem of 



plasma acceleration driven by radiation pressure in spectral lines provided that the flow is 
launched in the vicinity of a compact object. It was shown that it is important to include 
gravitational red-shifting in the calculations of the radiation force. The radiation pressure 
depends on d(j)/dr and dvjdr^ while in the case of 0-type star wind it depends only on dv/dr. 
Proximity to compact object, i.e. strong ionizing radiation, makes it difficult for heavy ions 
to survive against being completely stripped of electrons, so the efficiency of this additional 
mechanism in real situation depends on the simultaneous solution for the ionization balance 
and accounting on other possible mechanisms (such as clumping) to prevent over-ionization 
of the flow. 

We are also motivated by the current accumulating evidence for the gravitationally red- 
shifted narrow absorption features in many AGN spectra as well as by observations of gravi- 
tationally red-shifted absorption lines in the X-ray burst spectra of neutron stars. Potentially 
interesting example of the latter is presented by the observations of the gravitationally red- 
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shifted absorption lines in the X-ray burst spectra of EXO0748-676 (jCottam et aLl (J2002l )). 
If the interpretation is correct, these features are produced close to the compact object, 
where extreme conditions are coupled with high amplitude fluctuations of the radiation field 
and shortest dynamical time scale. As a result, these lines may be highly variable or/and 



transient. The non-detection of gravitation ally red-shi 



other good looking candidate GS, 1826-24 (JKong et al. 



t ed lin es in the observations of the 
(120071 )) may be an example. 



Given the diverse nature of these objects, it is desirable to consider a model with a 
minimum number of free parameters and in so doing we assume our outflow to be spherically- 
symmetric and exposed to the continuum radiation of a spherical core. 

Realistic models should consider a departure of the outflow and the radiation field from 
spherical symmetry. Bending of the photon's trajectories in the strong gravitational field of 
the compact object, may play some role. In our calculations we did not take into account 
because we were concerned with regions of the flow located approximately at radii > lOr^ . 

In our calculations the idea of the equal frequencies surfaces (EFS), plays a major role. 
The strong gravitational field changes the shape and locus of such a surface. Their topology 
is complex; the gravitational field strongly distorts the P-Cygni profile. Some branches of 
the EFSs in the red-shifted part of the spectra are found to be in front of the core, meaning 
the possibility for the absorption component to be observed as red-shifted, with respect to 
the emission. 

From numerical calculations, which are second order accurate in f /c, it is established 
that a superposition of Doppler and gravitational shifting of frequency can distort the P- 
Cygni profile in such a way that blue- and red-shifted absorption features are observed 
simultaneously. Often the red-shifted absorption line is superimposed on the emission wing. 
This effect is, of course, strongest close to the compact object. However, the emission part 
is also stronger there, and it may smear the absorption trough. Necessarily, the second 
absorption arises if the velocity profile is not very steep and at the same time the line forming 
region is situated within several tens of r^. However the former requirement is not crucial, if 
the the line opacity peaks close to the compact object. Profiles with more than one emission 
and two absorption features are possible. However, in a model of pure line scattering, those 
additional (i.e supplementary to W-shape profile) features would be considerably weaker. 

Different modifications of a "stellar" type velocity laws were adopted. Such velocity 
profiles describe stationary spherically-symmetric hydrodynamical flows. The well known 
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example of such is the v ~ Jl — 1/x profile. Note that different modifications of this law were 
obtained by different au thors using both theoretical and observational arguments (see e.g. 



Castor fc Lamerd (1l979l )). Thus we adopted a generalized form of this law, v ~ (1 — 1/x"^)'^ 



In case of non-magnetic accretion disk winds, the azimuthal component of the velocity 
quickly becomes much smaller than the radial one, and such a velocity profile describes the 
real distribution of the line-of-sight velocity reasonably well. 

During X-ray bursts, the situation can be more complicated. During the burst, the tem- 
perature at the stellar photosphere may quickly approach T ~ lO^ii' and the electron scat- 
tering occurs in the Klein-Nishina regime. The reduced cross section allows for higher ra- 
diation flux (locally below Eddington) to diffuse out to larger radii where temperature is 
smaller and the radiation flux is locally supper Eddington. A quasi-stationary wind may re- 
sult. Calculations of such winds must self-consistent ly account on the ra diation transfer and 



relatiy i stic corrections and demonstrate stellar type velocity profiles (JQuinn fc Paczynski 



(!l984J ): lNobihet all fll994J ).) 

To approximate the explosively accelerated plasma of X-ray bursters we adopt a linear 
velocity law. Making use of a Hu bble law, v = r/t d rastically simplifies the radiation transfer 



calculations in supernova shells (iKarp et al.l (119771 )). 

Two parameterized opacity laws and a linear, and generalized form of the velocity law, 
w ~ (1 — l/a;°^)°^ were considered. 

Our results for all considered velocity and opacity distributions show that in particular 
circumstances, i.e. proper velocity and opacity laws, strong gravity, the observed line profile 
consists of two absorption troughs separated by the emission component. The red-shifted 
absorption line can be weak and can be superimposed on the emission wing. Some of the 
strongest W-shape profiles were found for the linear velocity law. Surprisingly enough, no 



W-shape profiles are found for af~ wl — 1/x law. In this case the red-shifted absorption 
almost exactly compensates for the blue-ward emission peak, producing a sharp edge in the 
emission (Fig. [9]). 

However, in many other cases of different ai,a2, W-shape profiles are also found (Fig. [TO] 
- fTTl) . If the opacity is distributed as r ~ (1 — v/V°^Y, {k > 0), the red-shifted absorption 
line in most cases is superimposed on the red-shifted emission. However, if the opacity 
peaks at smaller velocities, r ~ l/v'^, distinct W-shape profiles are found again. Thus, the 
results suggest that the most favorable conditions realize when the opacity peaks at large 
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gravitational red-shifts and in such a case the velocity profile is probably less important. 
The further the distribution of the opacity is from that the more important is the velocity 
law and the relative importance of the Doppler effect. 

If some observed spectral feature is interpreted as being formed in the flow in a proxim- 
ity of a strongly gravitating compact object, then our results suggest that the assumption 
of the distribution of the opacity (optical depth law) is the most critical one. Such distri- 
bution should be calculated simultaneously with the ionization balance calculations in the 
background of the pre-calculated hydrodynamical model. 

The separation between red- and blue-shifted absorption features is a function only of 
the dynamics and relative importance of gravity in the line forming region. 

The most robust prediction of our model is the possibility of the double absorption 
trough as a result of a concurrence between gravitational red-shifting and Doppler effects. 
Our calculations suggest that the particular shape and intensity of the emission component 
which separates these absorption lines depend sensitively on assumed parameters. 

Gravitationally red-shifted absorption lines form in places close to the compact object. 
The emission component is necessarily formed in plasma occupying much larger volume. 
Thus, these features are formed in places which are possibly strongly separated in space. 
Thus, from the perspective of future observations, it would be interesting to look for cor- 
related variability of different components of the proflle between each other and with the 
continuum. 
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